Constraint on nuclear symmetry energy from nuclear charge 
dependent neutron skin thickness of nuclei 
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Abstract 

An alternative way to constrain the density dependence of the symmetry energy with the slope 
of the neutron skin thickness of nuclei which shows a linear relation to both the isospin asymmetry 
and the nuclear charge with a form of Z 2//3 is proposed. The linear dependence of the neutron 
skin thickness on the nuclear charge and isospin asymmetry is systematically studied with the data 
from antiprotonic atom measurement experimentally and the extended Thomas-Fermi approach 
incorporated the Skyrme energy density functional theoretically. An obviously linear relationship 
between the slope parameter L of the symmetry energy and the slope of the neutron skin thickness 
on the isospin asymmetry can be found by adopting 66 Skyrme interactions in the calculation. 
Combining the available experimental data, the constraint of 16 < L < 66 MeV on the slope 
parameter of the symmetry energy is obtained. The Skyrme interactions satisfying the constraint 
are selected. 
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The nuclear symmetry energy S(p) is the difference in energy per nucleon between pure 
neutron matter and symmetric nuclear matter, which is the key ingredient of the nuclear 
equation of state(EoS) for asymmetric nuclear matter. It governs the important properties 
of nuclei and neutron star. It also plays a significant role in the nuclear reaction dynamics 
and the stability of the phases within the neutron star and the interior cooling process 
within it. It is well known that the density dependence of S(p) for cold nuclear matter 
predicted from different models is extremely variant. Acquiring more accurate knowledge 
of the density dependence of the symmetry energy has become one of the main goals in 
nuclear physics at present and in the near future and has stimulated many theoretical and 
experimental studies lH5] . To characterize the density dependence of the symmetry energy, 
S(p) is expanded near saturation density (p ) as 
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The calibration of neutron skin thickness of nuclei defined by AR np = (r^) 2 — (f^) 2 has 
attracted a lot of attention in recent years because of the sensitivity of AR np to the density 
dependence of the symmetry energy. The calculations in either non-relativistic or relativistic 
mean-field models show a well-defined linear correlation between the AR np of heavy nuclei 
and the slope parameter L of the symmetry energy at the saturation density p-l8|. Thus, 
AR np of nuclei can be used as a powerful observable to constrain the density dependence of 
the symmetry energy at po and lower densities. The difficulty in the calibration of neutron 
skin thickness of nuclei stems from the difficulty in the measurement of neutron distribution. 
The main methods to measure the neutron distribution or neutron skin thickness include 



hadron scattering 



giant dipole 



12], 7r elastic scattering 



18M2C| and spin- dipole resonances 



131], antiprotonic atoms 14-12], excitation of 
221 ] on inelastic alpha scattering. Unfortu- 
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nately, the obtained values of AR np from different experimental method depend on the used 
analysis model and sometimes are not totally consistent with each other. It is also hard to 
judge the model dependence of the systematic error in different experimental method. The 



parity-violating electron scattering 
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24j will be a hopeful option to measure the neutron 



distribution with unprecedented precise of 1% in a model independent way. However, it is 
still not available up to now. In this case, it will be difficult to accurately and consistently 
constrain the symmetry energy by directly using the data of neutron skin thickness. We 
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notice that the AR np of 26 stable nuclei all over the periodic table 
have been accumulated from antiprotonic atom measurement. In Ref. 
of AR np on the isospin asymmetry 5 = (N — Z)/A for these 26 nuclei was extracted from 
experimental data of antiprotonic atom measurement, which reads 

AR np = (-0.03 ± 0.02) + (0.90 ± 0.15)5. (2) 

Recently, Warda et al. represented this relationship in droplet model with surface width 



dependence [251] . However it is already known that the antiprotons are only sensitive to the 
tail of the neutron distribution. An assumed shape for the neutron density is needed to 
extract the rms radius. Therefore the uncertainty in the value of AR np is unavoidable in 
this approach. 

In this letter, we suggest an alternative way to constrain the density dependence of the 
symmetry energy by the slope of a nuclear charge dependent linear relation of AR np vs 5 
based on the 26 experimental data from the antiprotonic atoms obtained up to now. In 
this way, only the tendency of AR np changing with the isospin asymmetry 5 is required 
to get the information of the symmetry energy. The systematic uncertainty due to the 
experimental method itself is therefore expected to be largely reduced. The more data of 
AR np for different nuclei with the same experimental method, the more accurate constraint 
on the density dependence of the symmetry energy can be obtained. One should take into 
account the nuclear charge dependent linear relation of AR np vs 5. Because the neutron 
skin thickness is the difference of the neutron and proton rms radii, it should depend not 
only on the symmetry energy but also on the Coulomb interaction which is far from clear. 

Let us first study the systematic behavior of the nuclear charge dependence of AR np . The 
approach of semi-classical extended Thomas- Fermi approximation H] up to the second 
order (ETF2) incorporated a potential energy density functional including the standard 
Skyrme energy density and Coulomb energy density with the Coulomb exchange term is 
applied in the calculations, in which the nuclear surface diffuseness is self-consistently taken 



into account. We calculate the proton and neutron density distributions of nuclei (See [28| for 



details) with a spherical symmetric Fermi functions by means of restricted density variational 



method 
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29h31|. With the density distributions determined in this way, the ground state 



properties such as the energy and the nuclear charge radii of series of nuclei have been 



calculated. The corresponding experimental data can be reasonably well reproduced 
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FIG. 1: Linear correlation between the neutron skin thickness and the isospin asymmetry for Ca, 
Ni, Zr, Sn, Yb, Pb and U elements. Scatter symbols denote the calculation results. Solid lines 
denote the linear fitting results with a form CzZ 2 ! 3 + C$5. 

Based on the calculated rms radii of proton and neutron, we get the neutron skin thickness 
AR np . The effective Skyrme interaction SLy4 is adopted in this work since SLy4 is very 



successful in describing the bulk properties and surface properties of nuclei [33|. 

With this approach we calculate the AR np for the isotopes with charge number in the 
range of 20< Z <92, only the even-even nuclei are taken into account. Fig.l shows the 
correlations between AR np and S for Ca, Ni, Zr, Sn, Yb, Pb and U isotopes. It is seen from 
the figure that the AR np and the isospin asymmetry 5 of nuclei are linearly correlated for 
Ca, Ni, Zr, Sn, Yb, Pb and U isotopes, respectively. The fitting lines are nearly parallel 
with each other and the slopes of AR np as a function of 5 for all selected elements is about 
1.10 fm per unit 8 within the range of < 5 < 0.24. The intercepts of the lines depends 
on the nuclear charge and roughly has a form of oc Z 2 ! 3 . We have varied the form of the 
charge number dependence from Z 1 ^ 3 to Z A I 3 and find that only the form of oc Z 2 ! 3 can 
describe the charge dependence of AR np . To understand the oc Z 2 ! 3 dependence of AR np , 
we investigate the relation between the squared charge mean radii and the charge number 
Z for the nuclei with same isospin asymmetry. We find that the squared charge mean radii 
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FIG. 2: Measured charge mean square radii (r 2 h ) 33] as a function of Z 2 / 3 for selected nuclei 
with the same 5. 

for the nuclei with same isospin asymmetry relate to Z 2 / 3 linearly. Fig.2 shows a monotonic 
increasing tendency of the experimental data of (r 2 h ) 



with Z 2 / 3 for nuclei with 5 — |, 
|, and 0.2, respectively. Obviously, increasing the charge number for the nuclei with same 
isospin asymmetry enlarges the (r 2 h ) and thus reduce the neutron skin thickness. 

According to the above investigation, we propose an empirical expression for the descrip- 
tion of the neutron skin thickness of a nucleus as 

AR np = C Z Z 2 ' 3 + C s d. (3) 

We have performed a \ 2 analysis in order to obtain the optimized C z and C$ from the 
experimental data of AR np from antiprotonic atom measurement. But the accumulated 
data of AR np are not enough to get a precise constraint on the parameters. Thus we obtain 
the parameters Cz and C$ by fitting the expression (2) obtained from experimental data of 
antiprotonic atom measurement. Firstly we use the data AR np for 40 Ca to determine the 
range of Cz = (—0.0041 ±0.0027)fm because for 40 Ca the 8 is equal to zero. Then we adopt 
Cz = —0.0041 fm to fit AR np for all 26 experimental data available up to now, in which 
the weight from error bar are taken into account to obtain the optimized C$- The fitting 




FIG. 3: (Color online) (a) Comparison of the results with (-0.0041 ± 0.0027)Z 2 / 3 + (1.05 ± 0.08)5 



14j. Results of DM are also 



to the experimental neutron skins from antiprotonic measurements 
shown. The dashed gray lines denote the upper and lower limits of Eq.(2) and the solid red lines 
denote those of Eq.(4). (b) to (e) Comparison of the results with the fitting linear formula to the 
experimental data from antiprotonic measurements for Sn, Zr, Te and Cd isotope chains. 

Cg of 1.05 fm is obtained. The range of C$ = (1.05 ± 0.08) fm is determined through two 
limits in the expression (2) extracted from experimental data. We finally obtain the neutron 
skin thickness as a function of both the nuclear charge and the isospin asymmetry of nuclei, 
extracted from the available data of AR np given by antiprotonic atom measurement as 

AR np = (-0.0041 ± 0.0027)Z 2/3 + (1.05 ± 0.08)5. (4) 

The uncertainty of C$ determined in this approach is considerably reduced from 0.15 fm in 
Eq.(2) to 0.08 fm when the charge dependence of neutron skin thickness is taken into account. 
The comparison between the AR np obtained with this expression and experimental data is 
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FIG. 4: (Color online) The correlation between L and the parameters Cz and Cs, respectively. 
The full circles are the calculated Cz and Cs for each Skyrme interaction. The solid black vertical 
lines denote the extracted value of Cz and Cs from experiment. The dashed black vertical lines 
denote the extracted upper and lower limits of them. 



given in Fig. 3 (a). The results of droplet model (DM) 25| are also shown in Fig. 3 (a) for 
comparison. One can see that Eq. (4) can reproduce most of the experimental data well and 
also is in good agreement with the DM calculations. In the sub-Figs (b),(c),(d),(e) of Fig.3, 
we pick out Sn, Zr, Te, Cd isotope chains from the available 26 experimental neutron skin 
thickness data and compare with the results of Eq.(4). One can see that the experimental 
data are reproduced reasonably well. 

Now let us come to explore the density dependence of the nuclear symmetry energy from 
the neutron skin thickness of nuclei AR np correlated with both charge and 6 of nuclei. We 
perform a systematic calculations within the ETF2 approach by using 66 sets of Skyrme 
interactions like SLy series, SkT series, v series, Skz0~4, MSkl~6, SkSCl~4, SkM, SkM*, 
SkMl, SkMP, SII ~ SVII, SUP, SGI, SGII, Zs, Es, Gs, Rs, FitB, BSkl, RATP, SKRA. The 
slope parameter L of the bulk symmetry energy covers a widely range from nearly —30 MeV 
to 100 MeV for all the 66 Skyrme interactions. For each Skyrme interaction, we calculate 
the neutron skin thickness of the nuclei for which the neutron skin thickness data from 
antiprotonic atom measurement are available. Then we get the parameters Cz and Cs in 
Eq.(3) for each Skyrme interaction. The relationships between the slope parameter L and 
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Gz y Cs thus can be obtained. Fig. 4 displays the correlations between the slope parameter 
L of the symmetry energy and the parameters Cz (left panel) and Cs (right panel). The 
full circles are the calculated Cz and Cs for each Skyrme interaction. From the left panel 
we can see that the parameter Cz changes from —0.006 to —0.003 fm which is totally 
in the range of [—0.0068, —0.0014] fm extracted from the experimental data, and shows 
an anti-correlation with L of Skyrme interactions. In the right panel of Fig.4, a clearly 
linear increasing correlation between L and Cs can be observed. The slope parameters 
L for different Skyrme interactions varies with Cs in an area limited by two lines L = 
(-117.93 ± 2a) + 151.21(7,5 (a = 6.54) which are plotted in Fig.4(b) with dashed red lines. 
With the range of Cs = (1.05 ± 0.08)fm extracted from the experimental data, the slope 
parameter L is constrained in 16 < L < 66 MeV. Within the 66 Skyrme interactions 
we adopted, the Skyrme interactions satisfying this relation are SLy series, SkM, SkM*, 
MSkl~2, Skz0~l, SkT3, SkT6~9, SII, SUP, SGII, RATP. The corresponding symmetry 
energy coefficients at the saturation density po are from 28 MeV to 32 MeV for these selected 
Skyrme interactions. The uncertainty of L obtained from the extracted range of Cs is 
much smaller than that from Cz- It means that the correlation between AR np and isospin 
asymmetry 5 of nuclei is more sensitive to the density dependence of the symmetry energy 
than that between AR np and nuclear charge. The obtained constraint on the slope parameter 
L of the symmetry energy is in consistent with that obtained from the double neutron-proton 
ratios and isospin diffusion in heavy ion collisions [lj]. 

In summary, we propose an alternative way for constraining the density dependence of 
the symmetry energy by means of the relation of the neutron skin thickness to the nuclear 
charge and isospin asymmetry of nuclei. We show that the neutron skin thickness depends 
on the isospin asymmetry and the charge of nuclei with the form of Z 2 ^ 3 . By fitting the 
available data of neutron skin thickness for a series of nuclei obtained from antiprotonic atom 
measurement, the parameters for the correlation between AR np and the nuclear charge (C z ) 
and isospin asymmetry (Cs) can be extracted. Within the framework of the semi-classical 
extended Thomas-Fermi approximation together with the Skyrme energy density functional 
and Coulomb energy density, we systematically calculate the neutron skin thickness with 66 
Skyrme interactions for the 26 nuclei, for which the experimental data are available, and get 
the parameters C z and Cs for each Skyrme interaction, respectively. Based on the clearly 
linear correlation between the slope parameter of the symmetry energy and the slope pa- 
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rameter C$ of the neutron skin thickness, we obtained the constraint on the slope parameter 
of the symmetry energy, i.e. 16 < L < 66 MeV, by using the extracted parameter C$ from 
the experimental data. The Skyrme interactions satisfying this constraint are selected. It 
should be stressed that the spherical symmetry Fermi distribution for proton and neutron 
density in our calculation is in accordance with that in the analysis in the antiprotonic mea- 
surement. It is more suitable and consistent to extract the parameters C z and C$ with the 
data from antiprotonic measurement. We suggest more neutron skin thickness for different 
element with the same isospin asymmetry to be measured experimentally, so as to extract 
a more accurate dependence of the neutron skin thickness on the nuclear charge and the 
isospin asymmetry. With the increasing of the number and the accuracy of data for neutron 
skin thickness of nuclei based on the antiprotonic atom measurement in future, we believe 
more accurate constraint on the density dependence of the symmetry energy can be obtained 
with our proposed method. 
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